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We investigate some exact static and dynamic properties of one-dimensional fermionic Tonks- 
Girardeau gases in tight de Broglie waveguides with attractive p-wave interactions induced by a 
Feshbach resonance. A closed form solution for the one-body density matrix for harmonic trapping 
is analyzed in terms of its natural orbitals, with the surprising result that for odd, but not for even, 
numbers of fermions the maximally occupied natural orbital coincides with the ground harmonic 
oscillator orbital and has the maximally allowed fermionic occupancy of unity. The exact dynamics 
of the trapped gas following turnoff of the p-wave interactions are explored. 
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The spin-statistics theorem, according to which iden- 
tical particles with integer spin are bosons whereas those 
with half-integer spin are fermions, breaks down if the 
particles are confined to one or two dimensions. Real- 
ization of this fact had its origin when it was shown by 
one of us that the many-body problem of hard-sphere 
bosons in one dimension can be mapped exactly onto that 
of an ideal Fermi gas via a Fermi-Bose (FB) mapping. It 
is now known that this "Fermi-Bose duality" is a very 
general property of identical particles in one-dimension 
(ID). In recent years this subject has become highly rel- 
evant through experiments on ultracold atomic vapors in 
atom waveguides 0, |1 0, i, |1 i, i, E [H E [S [11 . 
Exploration of these systems is facilitated by tunability 
of their interactions by external magnetic fields via Fes- 
hbach resonances jl5l |. In fermionic atoms in the same 
spin state, s-wave scattering is forbidden by the exclusion 
principle and p-wave interactions are usually negligible. 
However, they can be greatly enhanced by Feshbach reso- 
nances, which have recently been observed in an ultracold 
atomic vapor of spin-polarized fermions Additional 
resonances are induced by tight transverse confinement in 
an atom waveguide. Of particular interest is the regime 
of low temperatures and densities where transverse os- 
cillator modes are frozen and the dynamics is described 
by an effective ID Hamiltonian with zero-range mterac- 
tions FitI 11^. a regime already reached experimentally 
[lIlIMinfElIIllli- Transverse modes are still virtu- 
ally excited during collisions, leading to renormalization 
of the effective ID coupling constant via a confinement- 
induced resonance. This was first shown for bosons 
and recently explained in terms of Feshbach resonances 
associated with bound states in closed, virtually excited 
transverse oscillation channels , and recently Granger 
and Blume 20] have shown that the same phenomenon 
occurs in spin-polarized fermionic vapors. 

For a system of bosons with hard core repulsive inter- 
actions, a Tonks-Girardeau (TG) gas, the FB mapping 
relates the ground state for the system to the ground 
state of a free Fermi gas [i^. We have recently pointed 
out 13] that in the case of the spin-aligned Fermi 



gas, the generalized FB mapping USUI Hi can be ex- 
ploited in the opposite direction to map the "fermionic 
TG gas" 113,13]' ^ spin-aligned Fermi gas with strong ID 
atom-atom interactions mediated by a 3D p-wave Fesh- 
bach resonance, to the ground state of the trapped ideal 
Bose gas. In this Letter we study some new and exact 
static and dynamic properties of fermionic TG gases. 

Spin-aligned Fermi gas: A magnetically trapped, spin- 
aligned atomic vapor of spin-^ fermionic atoms in a tight 
waveguide is magnetically frozen in the spin configuration 
Ti ' ' ■ Ta^j so the space-spin wave function must be spa- 
tially antisymmetric, s-wave scattering is forbidden, and 
the leading interaction effects at low energies are deter- 
mined by the 3D p-wave scattering amplitude. For a ID 
gas of free spin-aligned fermions with harmonic trapping 
of frequency u the Hamiltonian is 
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and the A'^-fermion ground state V-Fola^i, • ' ' ^xn) is a 
Slater-determinant of the A^ lowest HO orbitals (f>n{x) = 
e-Q'/2ff„(Q)/7ri/4a;i/2y2^ ^ ^ith iJ„(g) th e Her- 
mite polynomials and Q — x/xosc where Xosc — y^h/muj 
is the longitudinal oscillator width. The Fermi-Bose map- 
ping method 1] maps this spin-aligned A'^-fermion ground 
state to the A^-boson ground state of a system of ID im- 
penetrable bosons (TG gas) leading to an explicit expres- 
sion for the A^-boson ground state |2^ . 

Here we consider the opposite limit of the "fermionic 
TG gas" 113,01 1 * spin-aligned Fermi gas with strong ID 
atom-atom interactions induced by a 3D p-wave Feshbach 
resonance |0. Granger and Blume derived an effective 
one-dimensional K-matrix for longitudinal scattering of 
two spin-aligned fermions confined in a single-mode har- 
monic atom waveguide (2(]|. It can be shown [13,111111 
that in the low-energy domain the K-matrix can be re- 
produced, with a relative error as small as C(fcf ), by the 
contact condition 



i^pixje = 0+) = -tpF{xji = 0-) = -a'[jjil)p{x.ji. 



= 0±) 
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where the prime denotes differentiation, 

«fz5 = ^[i + i2(v4)ic(-iAi)r^ (3) 

is the one-dimensional scattering length, a± = y^h/ ^u!± 
is the transverse oscillator length, Vp = — 
— limfc^o tan(5p(fc)/fc^ is the p-wave "scattering volume" 
[isf . Qp is the p-wave scattering length, ^(—1/2,1) = 
-C(3/2)/47r = -0.2079. . . is the Hurwitz zcta function 
evaluated at (—1/2, 1) 26], and /i = m/2 is the reduced 
mass. The expression has a resonance at a nega- 
tive critical value Vp"* /a^ = —0.4009 • • • . In accordance 
with 1(21, the low-energy fermionic wavefunctions are dis- 
continuous at contact, but left and right limits of their 
derivatives coincide, and following Ref. j2ll | we assume 
the same here. For an interatomic interaction of short 
but nonzero range xq, the fermionic relative wave func- 
tion vanishes at xji — 0; Eq. Q refers to the xq 0+ 
limit of the exterior wave function just outside the range 
of the interaction H^, H^. 

The contact condition (01 is generated by the 
ID interaction pseudopotential operator |3 
vfnt = 9iDlli<j<ii<N^' i^3f)9jt where = 
{\/2)[dxjik\xj=xi+ - 9x^'0Uj=x^^] with_effe£tive ID cou- 
pling constant gfjj = h^afjj/^i [la, IH llg. This can 
be compared with the bosonic ID coupling constant 
gfjj = —ti^/iiafjj llSHi- The dimensionless fermionic 
coupling constant is = —'mgiD'n/fi^. Note that the 
density n is in the numerator, whereas it is in the de- 
nominator of the bosonic analog 7b ~ mgfjj/nh^. 

Fermionic TG gas: The fermionic TG gas regime 
is 7f S> 1 valid in the neighborhood of the reso- 
nance where the denominator of Eq. vanishes. 
The ID scattering l eng th is invariant under the FB 
mapping P, HO, '0b — A{xi,- ■ ■ ,xn)4'f with 

A{xi,- ■ ■ ,xn) = I\i<j<e<N^S^i^e ~ ^j) ^^^^ "unit an- 
tisymmetric function" employed in the original discovery 
of fermionization 1], so we shall henceforth indicate it 
by fliD — afjj = afjj. Since Vp must be negative to 
achieve a resonance, the fermionic TG regime with at- 
tractive p-wave interactions auj — > —oo and — > —oo 
is reached as Vp approaches the resonance through neg- 
ative values, implying an interaction-free exterior wave 
function. More generally, the spin-aligned Fermi gas 
maps to the Lieb-Liniger Bose gas |2^ with fermionic 
and bosonic coupling constants inversely related jl^ ac- 
cording to g^jy = —h'^/fi'^g[j~, > 0. The corresponding 
dimensionless coupling constants 75 = mgf^Jrifi^ > 
and 7f = —mgij^n/fi? > satisfy ^bIf — 4 !22|. 

Natural orbitals: The A^-atom ground state 
of the longitudinally trapped fermionic TG gas 
maps to the trapped ideal Bose gas ground state 
(7b — 0) for which all N atoms are Bose-Einstein 
condensed into the lowest HO orbital 4>q{x): 

Ippixi,- ■ ■ ,XN;t = 0) = A{xi,- ■ ■ ,XN)l\j = i(f>Q{Xj). 



Since = 1, all properties expressible in terms of 
including the density profile, are the same as 
those of the trapped Bose gas ground state Ylj^i 4>o{xj)- 
However, properties not so expressible, such as the 
momentum distribution, differ dramatically, as in 
the case of the usual bosonic TG gas. The ap- 
propriate tool for studying such differences is the 
reduced one-body density matrix (OBDM) pi{x,x';t) = 

N J ipF{x,X2, ■ ■ ■ ,XN;t)^p*p{x' ,X2, ■ ■ ■ ,XN;t)dx2 ' ' -dxN- 

The normalized natural orbitals satisfy XjUj{x) = 
/_ dx' pi{x,x';t = 0)uj{x'), with eigenfunctions Uj{x), 
j = 0,1,2,..., and eigenvalues or orbital occupations 

< Aj < 1. For a spatially uniform bosonic TG 
gas of length L the natural orbitals are plane waves 
"i(^) — e^"^""'^^!^ I^Tj and the distribution of occupa- 
tions Aj reflects a momentum distribution function with 
low-momentum peaking and broad wings, whereas for 
a free Fermi gas the ground state is a filled Fermi sea 
with \j = Lj_2 = 0^, . . . ,N — 1, all other occupations 
being zero 0, Q, I23 . Although natural orbitals are 
much-used in quantum chemistry, the first application 
to degenerate gases is due to Dubois and Glyde js^l . 

In the case of the longitudinally trapped 
fermionic TG gas, the OBDM can be eval- 
uated in closed form for all N 31] yielding 
piix,x';t = 0) = N(l)o{x)(l)oix')[F{x,x')]^'^ where 
F{x,x') - j:^^sgn{x - x")sgn{x' - x")^l{x")dx" = 

1 — ]erf(Q') — erf((5)], with erf((5) the error function |3^ . 
We have numerically diagonalized the OBDM to obtain 
the natural orbitals and occupations for a variety of N. 
In Figs.^a) and (b) we show the occupations Xj versus 
orbital number j for N — 2 and TV = 3, respectively, 
as representative of our findings for even and odd 
numbers of atoms. (The solid line between data points 
is included as an aid to the eye). For TV = 2 we see that 
the occupancy distribution shows a staircase structure 
where X2q — A2g+i, g = 0, 1, 2 . . . , the occupancies being 
monotonically decreasing Xq > Ag+i. Furthermore, the 
occupancy distribution extends well beyond the Fermi 
sea for the corresponding free Fermi gas shown as the 
dashed line, reflecting the strongly interacting nature of 
the fermionic TG gas . For = 3 shown in Fig. ^b) 
a similar staircase structure is seen with the exception 
that Ao = 1, that is the natural orbital uo{x) has the 
maximal occupancy allowed by the exclusion principle, 
and X2q+i ~ X2q+2, g = 0, 1, 2 . . . . 

The natural orbitals are structurally similar to the HO 
orbitals and share the property that the orbital number 
j is equal to the number of field nodes. To examine the 
properties of the maximally occupied natural orbital we 
used a Rayleigh-Ritz variational procedure with a Gaus- 
sian trial solution uo{x) = exp(— a;^/2w^)/7r-^/'*-\/w to ob- 
tain the corresponding occupancy Ao (w) numerically. We 
have found that for odd N the maximal occupancy ob- 
tained by varying w is Xo{w = Xosc) = 1 for w = Xosc, 
so that the maximally occupied natural orbital uo{x) co- 
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(a) 




Natural orbital number j 

FIG. 1: Occupation distribution \j versus natural orbital 
number j for a) A*' = 2, and b) A*' = 3. In each case the 
dashed line shows the corresponding result for a free Fermi 
gas. 



incides with the ground state (/)o {x) of the single-particle 
HO. We have verified this numerical result by proving 
analytically, by tedious rearrangement and summation 
of series, that Aq = 1 assuming uo(x) = 0o(a;) for odd 
N , an exact and surprising result. In contrast, for even 
iV the maximum occupancy is obtained numerically for 
w = Xosc/'^i the occupancy increasing with N but re- 
maining less that unity, and comparison of the numerical 
results and Gaussian solution shows that it is an approxi- 
mation to the maximally occupied natural orbital, albeit 
a very good one. 

Dynamics following turnoff of interactions: We next 
examine the exact dynamics of the fermionic TG gas 
that ensue when the interactions are suddenly turned 
off at t = by detuning the external magnetic field. 
For t < the OBDM is given by that for the fermionic 
TG gas pi{x,x';t < 0) = ^-\jUj{x)uj{x'). For later 
times, when the interactions are turned off, the evo- 
lution of the OBDM is given by ih{d/dt)pi{x,x';t) — 
[H{x) ~ H{x')]pi{x,x']t) js^, where H{x) is the single- 
particle HO Hamiltonian. Even though this equation has 
no interactions the subsequent dynamics are a direct con- 
sequence of the strong many-body correlations present in 
the initial OBDM. Formally, the solution for the OBDM 
can be written in terms of the HO orbitals that serve as 
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FIG. 2: Scaled density profile Xoscp{Q,t = T/2) after half a 
period for a) A'^ = 2 and N = 8 (dashed line), and h) N = 3 
and N — 9 (dashed line). 

natural orbitals for the non-interacting gas 
pi{x, x';t>0)=J2 C™„0„(x)0„(x')e-'("-'")'^*, (4) 

where Cnm = J2j I dxun{x)4>j{x) J dxum{x)(t)j{x), 
and the symmetry of the orbitals dictates that Cnm is 
non-zero only for n and m both even or both odd. Thus, 
the OBDM evolution following the turnoff of the interac- 
tions is periodic in time with period T — t:/uj. We have 
solved the equation of motion for the OBDM numerically 
to obtain the exact dynamics of the initial fermionic TG 
gas, and here we present results for the evolution of the 
density profile p{x, t) = pi{x, x; t). 

Figure [21 a) shows the scaled density profile 
XoscPiQit — T/2) after half a period, when the 
on-axis density reaches its lowest value, for = 2 and 
N = 8 (dashed line). In each case we see that the density 
expands and the center {Q = 0) develops a minimum. 
This may be intuited as follows: Before the interactions 
are turned off the width of the density profile is ~ Xosc, 
narrower than the width for a trapped free Fermi gas 
^ XoscVN ■ This difference in widths arises from the fact 
that the Fermi degeneracy pressure that sets the width 
of the free gas in conjunction with the HO trapping po- 
tential is countered by the attractive p-wave interactions 
in the fermionic TG gas. However, when the interactions 
are turned off the Fermi degeneracy pressure, which 
increases with N , acts to broaden the density profile. 
For times t > T/2 the density profile begins to contract 
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again under the action of the harmonic single particle 
potential, the initial condition being recovered at i = T. 

Figure [2Ib) shows the corresponding results for = 3 
and = 9 (dashed line), and similar broadening is seen 
at i = T/2. In this case, however, a density maximum 
is maintained at the center of the trap. The persistence 
of an on-axis density peak for odd N may be understood 
as follows: For odd TV the particles are symmetrically 
placed around the center of the trap with one particle 
at the center. When the interactions are turned off the 
particles displaced from the center will move outwards in 
reaction to the Fermi degeneracy pressure. However, the 
center particle, which is absent for even N , will be pinned 
at the center by the equal but opposite Fermi degeneracy 
pressures acting on it due to the displaced particles, and 
this maintains the on-axis density peak for odd N . The 
persistence of a central peak may also be related to the 
fact that for odd N the maximally occupied natural or- 
bital uq{x) for the fermionic TG gas with interactions is 
identical to the HO orbital 0o {x) which serves as a natu- 
ral orbital with no interactions. Thus, the scaled density 
profile associated with this maximally occupied orbital 
is given by XoscPoiQ) = CooXosc4'l{Q) = exp{-Q'^)/y/TT, 
which has a maximum on-axis, integrates to unity and 
represents the centrally pinned particle, and persists af- 
ter the interactions are turned off. For the = 3 
simulation in Fig. |2Ib) Xoscp{0,T/2) ~ 0.6, whereas 
XoscPoi^) — 0-56, showing that the central peak is largely 
accounted for by the maximally occupied natural orbital. 
For larger odd N the central maximum has a larger con- 
tribution from higher natural orbitals but the peak still 
persists; see the dashed line in Fig. I^Jb). 

In summary, the OBDM of a harmonically trapped 
fermionic TG gas has been analyzed in terms of its natu- 
ral orbitals and occupations. For an odd number of atoms 
the maximally occupied natural orbital was found to co- 
incide with the single-particle HO ground state. Due to 
the rarity of exact results for strongly interacting systems 
it is important to suggest observable consequences. We 
have explored the exact dynamics following turnoff of the 
p-wave interactions, finding that for odd N the on-axis 
density remains a maximum as the cloud expands due to 
Fermi degeneracy pressure, in contrast to even N where 
the density develops a minimum. 
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